Vibration Spectra of the m-Tree Fractal 

Daniel J. Ford and Benjamin Steinhurst* 
December 15, 2008 



Contact: 
steinhurst @ mat h . uconn.edu 
Department of Mathematics 
University of Connecticut 
Storrs CT 06269 USA 

Abstract 

We introduce a family of post-critically finite fractal trees indexed by 
the number of branches they possess. Then we produce a Laplacian opera- 
tor on graph approximations to these fractals and use spectral decimation 
to describe the spectrum of the Laplacian on these trees. Lastly we con- 
sider the behavior of the spectrum as the number of branches increases. 
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1 Introduction 

There have been many studies using spectral decimation to calculate the spec- 
trum of a Laplacian operator on finitely ramified fractals [TJ [5J [JJ] [TBI EZ] • These 
authors have studied many fractals in both specific cases, in families, or in gen- 
eral. We shall consider a Laplacian for a family of fractals, the m-Branch Trees 
and in one general calculation find the spectrum of Laplacians on the entire fam- 
ily. This work is done via the spectral decimation process by which the spectra 
of Laplacians on graph approximations are used to calculate the spectrum of a 
limiting Laplacian on the fractal or on an infinite graph. Our notation is slightly 
complicated by the fact that we do the spectral decimation calculations for an 
arbitrary m, so effectively doing every fractal in this family at one time. Trees 
have been a well studied topic in fractal literature, for example Vicsek sets in 
PHOT] and Dendrites in 0. 

In [10], Kigami set out the frame work by which post-critically finite frac- 
tals can be thought of as abstract spaces independent of any ambient space. 
This is accomplished by labeling a point by an "address" determined by the 
cell structure of the fractal. The space of addresses is the fractal. Because of 
this point of view it is natural to consider graph approximations to the fractal 
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as truncations of these addresses. Kigami also proved that the graph Lapla- 
cian operators on the approximating graphs converge to a Laplacian on the 
full fractal. Strichartz, [15], builds upon Kigami's point of view and considers 
differential equations on post-critically finite fractals using the types of opera- 
tors that Kigami constructed. Strichartz included a discussion of the spectral 
decimation method for calculating the spectrum of these Laplacians from the 
spectra of the graph Laplacians. 

The m-Branch Trees, m 3, are post-critically finite fractals whose approx- 
imating graphs are defined for m 3 and are constructed with m-simplices. We 
run through the calculations of the spectrum for the 3-branch tree in Section [¥] 
before doing the general case for the sake of having concrete matrices to work 
with. Then in Section [5] we do the general case and finally in Section [6] we 
observe what the behavior of the spectra of the Laplacians on the m-branch 
trees is as the number of branches increase. But first we describe the trees in 
Section [2] and the spectral decimation method in Section [3l 
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2 m-Branch Tree Fractals 

The m-branch tree fractal, F mi is a post-critically finite fractal with m defining 
contraction mappings. The zero-level graph approximation, V^Oj consists of a 
complete graph of m vertices, when m = 3, V3.0 is a triangle, when m = 4, V^o 
is a tetrahedron. The iterated function system that generates the fractal scales, 
duplicates, and translates the simplex to m simplices sharing a common point at 
the epicenter of the previous simplex and with each vertex from V m fl as a vertex 
of one of the new simplices, this is the graph V m> i. This process is iterated and 
the countable set of vertices is completed in the effective resistance metric to 
form a tree with m branches. With the graphs V m fl and V mi i described we refer 
to the appendix of [7] for the detailed construction of the post-critically finite 
fractal and proofs of its properties. 

Below we construct a Laplacian operator, A m , on F m as a limit of Laplacians 
A m; „ on V m . n from which we will calculate the spectrum of A m by spectral 
decimation. All of the calculations that are done in the spectral decimation 
process do not depend on any particular embedding into a Euclidean space and 
in fact when we calculate the Hausdorff dimension of these space in Section [5] 
we will use the intrinsic effective resistance metric. This metric capitalizes on 
the analogy between graphs and electrical networks to calculate the effective 
resistance between two points as if the edges of the graph were resistors and the 
vertices nodes in an electrical network. This analogy has been used by many 
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authors but a useful reference for the mechanics of these calculations and a more 
extensive bibliography is [15]. It is a straight forward calculation to show that 
in the effective resistance metric the contraction mappings forming the trees 
have contraction factor one half for any to. 

There has recently been developed a connection between fractals and certain 
groups, most specifically cellular autonoma and iterated monodromy groups. 
Nekrashevych has written a survey |12j exploring self-similar groups such as it- 
erated monodromy groups involving Schreier graphs as pictoral representations 
of the groups. In [3] and [I] there is consideration of a particular self-similar 
group, which is a subgroup of the automorphism group of a rooted p— tree. This 
particular self-similar group is sometimes called the Gupta-Fabrykowski group. 
Bartholdi, Grigorchuk, and Nekrashevych explore the limits of the Schreier 
graphs as fractals combining algebraic and analytic perspectives. The Schreier 
graphs for the Gupta-Fabrykowski groups are the same as the V m ,n graphs for 
the TO-Branch Tree when m is a prime. 

Figures [JJ and [2] gives a visual guide to what the iterated function system 
does in the to = 3 case. Note that all triangles are equilateral in the intrinsic 
effective resistance geometry but when embedded in R 2 they would overlap if 
they were drawn as equilateral in the geometry of R 2 . Since the theoretical 
machinery that we use in this paper does not depend on a particular embedding 
into Euclidean space we won't be using the geometry of any R n so the issue of 
overlap is not relevant. 




V => Vl 

Figure 1: The progression of the three-branch tree from its V3,o to 1/3,1 network 

The Laplacian operator we use is analogous to the second difference op- 
erator on a continuous space because it takes the average value of the func- 
tion at neighboring points and subtracts that from the value of the function, 
A m . n f(x) = f(x) - Y,x~y dcg[>) • Where the sum is over neighboring points in 
the graph and deg(ir) is the number of neighbors that x has. Because the ma- 
trix representation of this operator on a graph has ones on the diagonal this 
is known as a probabilistic Laplacian. We construct a matrix representation 
for the Laplacian A„ i: „ associated to the graph V m , n as in [TJ [2j [T5l [16j HZ]- 
Let M„ l n denote this matrix which is the level- n graph approximation to the 
Laplacian on F m . The on-diagonal entries of M m>n are 1 as already noted. The 
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Figure 2: The graph V^, note that the branches only connect to each other at 
the center point. 



off-diagonal entries are: 

• if Xi and X j 0X6 not connected; 

• de ~( x ) if Xi and Xj are connected. 

Kigami proved that the sequence of operators p n M n converge in an appropriate 
sense to an operator A m on F m in [101 QT] where p is the "energy renormaliza- 
tion constant." If we view the distance between adjacent vertices always being 
constant as the level of approximation increases then we get an infinite graph 
instead of a compact fractal where p = 1 , if the graphs are scaled at each step so 
that the limiting fractal is bounded then p > 1. We will calculate the spectrum 
of the bounded operator limn-joo M m ^ n = A m and will not be concerned with 
p. Since there are deg{xi) points connected to Xi and the diagonal entry is 1, 
then we have the property that the row sums of M m ^ n are always equal to zero 
which is a very useful aid in the calculations. 



3 Spectral Decimation 

Spectral decimation is a means of extending the eigenvalues for the n-level 
Laplacian to those of the n + 1-level Laplacian. Used inductively, this means 
only the zero- and first-level Laplacians must be explicitly calculated. The 
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process, which dates back to [BKHUH], is presented in terms of the calculations 
shown in [2] and discussed in more abstract detail in [lj [T4J, [15] . We first write 
the level-one Laplacian as a block matrix. For the time being we drop the m 
from the notation for simplicity's sake. 

M = M t = ( % d\ (3 - 1} 

where A = I m since the boundary points are not neighbors in V\. The Schur 
Complement of the matrix M — Iz is the matrix-valued function: 

S(z) = (A-z)-B(D-z)- 1 C. (3.2) 

When v — (vq, v[) t is an eigenvector of M with corresponding eigenvalue z, the 
eigenvalue equation can be written as 



A B \ ( v \ ( v 
C D 



(3.3) 



This system of linear equations yields two equations: v[ = —(D — z) Cvo, 
provided that z (fi o~{D), and (A— z)vo+Bv[ = 0. It is worth noting that the map 
— {D — z)~ 1 C takes the eigenvector vq on the boundary points and determines 

the values required at the interior vertices so that [ ^ J is an eigenvector 



on V\. This is what actually does the eigenvector extension mentioned below. 
When the two equations are combined with the Schur complement they imply 
that 

S{z)v = 0. (3.4) 
If v is also an eigenvector of M with corresponding eigenvalue Zq, then 

(M - z )v a = 0. (3.5) 

By combining l|3.4p and (|3.5p . and allowing zq = R(z) for a rational function 
R(z), we have, 

S(z) = <f>(z) {M - R(z)) (3.6) 

where <j){z) and R{z) are scalar-valued rational functions whose existence is 
proved in [TJ [SJ [TJ1 [T5] . It can be observed from the proof of their existence that 
cj)(z) is dependent on 51.2 and R(z) is dependent on both Si t i and 4>{z). The 
role of R(z) is as an "eigenvalue projector" taking eigenvalues of one level and 
projecting them down to a lower level's eigenvalues. To calculate these functions 
we only need two entries from S(z). The formulae are: 

Si i 

4>{z) = —{in — \)S\ 2 and R{ z ) = 1 tt^t- (3-7) 

<P(z) 

These equations are a priori rational functions on C so the locations of 
possible zeroes and poles require special attention. If they are off the real 
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axis they won't be an issue since A m are bounded symmetric operators, only 
zeros and poles on the real axis are of concern. The location of these possible 
zeroes and poles are what we call exceptional values. The exceptional values 
are E(Mq,M) = a(D) U {z : <fi(z) — 0}. These values are exceptional because 
(D — z) is not invertible when z e cr(D), which causes (|3.4p to not be defined, 
and if <j>{z) = then either (|3.4p and/or (|3.5p fail to be defined, i.e. the poles 
of R and <j> and zeroes of 4>. 

It is the rational functions R(z) and 4>(z) that are the main tools used 
throughout this paper. Proposition 15.41 will produce a general form for the 
values of S^i and Si, 2 to obtain the functions <f>{z) (Corollary 15. 5[) and R{z) 
( Corollary in the general m-branch case. The function R{z) projects level- 
n + 1 eigenvalues to level-n eigenvalues, so the primary use of R(z) is to take 
inverse images under it. We will call the depth- n eigenvalue z the offspring of 
the depth-n + 1 eigenvalue R^ 1 (z). 

The process of extending eigenvalues from one approximation to the next 
is summarized in the following proposition from [T]. Above an extension map 
was mentioned, —(D — z)~ 1 C, from this extension map the eigenvectors are 
extended providing that z is not an exceptional value but is an eigenvalue at 
the next higher level. To find those eigenvalues at the next level we will use this 
proposition. Denote by mult^z) the multiplicity of z as an eigenvalue of D. 
Similarly, mult„(z) is the multiplicity of z as an eigenvalue of M n . Again the 
subscript m has been removed to simplify notation, so everything here is for a 
given value of m. 

Remark 3.1. Here and throughout, dim„ is the dimension of the function space 
on Vn, for a given m. Since V n is a finite collection of points, the dimension of 
the space of functions on V n is just the number of points in V n . 

Proposition 3.2. 

1. If z£ E(M Q ,M), then 

mult„(z) = mult„_ 1 (i?(z)), (3.8) 

and every corresponding eigenfunction at depth n is an extension of an 
eigenfunction at depth n — 1 . 

2. If z ^ o~(D), (f>(z) = and R(z) has a removable singularity at z, then 

mult„(z) = dim„_i, (3.9) 
and every corresponding eigenfunction at depth n is localized. 

3. If z £ o~(D), both 4>{z) and <p(z)R(z) have poles at z, R(z) has a removable 
singularity at z, and 4-R(z) 5^ 0, then 

mult„(z) = m n ~ rmilt£) (#) — dim„_i + mult n _i(i?(z)), (3.10) 

and every corresponding eigenfunction at depth n vanishes on V n -\. 
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4- If z G o~(D), but <j>(z) and <f)(z)R(z) do not have poles at z, and <j>(z) ^ 0, 
then 

mult„0) = m n-1 multr>(,s) + mult„_i(i?(z)). (3.11) 

In this case m n ~ 1 mult rj(z) linearly independent eigenf unctions are local- 
ized, and mult„_i {R{z)) more linearly independent eigenf unctions are ex- 
tensions of corresponding eigenfunction at depth n — 1 . 

5. If z G o~(D), but 4>{z) and <p(z)R(z) do not have poles at z, and 4>{z) = 0. 
then 

mult„(z) = r7i" _1 multz)(z) + mult„_i(E(z)) + dim„_i (3.12) 

provided R(z) has a removable singularity at z. In this case there are 
m™~ 1 mult£)(z) +dim n _i localized and mult„_i(-R(z)) non-localized corre- 
sponding eigenfunctions at depth n. 

6. If z G o~(D), both 4>(z) and 4>(z)R(z) have poles at z, R(z) has a removable 
singularity at z, and -j^R(z) — 0, then 

mult„(z) = mult n _i(i?(z)), (3.13) 

provided there are no corresponding eigenfunctions at depth n that vanish 
on V n -\. In general we have 

mult„(z) = m n ~ 1 mu\t D (z) - dim„_i + 2mult„_i(i?(z)). (3.14) 



7. If z ^ (t(D), 4>(z) = and R(z) has a pole z, then mult n (z) = and z is 
not an eigenvalue. 

8. If z G c(D), but 4>{z) and 4>(z)R(z) do not have poles at z, 4>{z) — 0, and 
R{z) has a pole z, then 

mult„(z) = m n - 1 nm\t D (z) (3.15) 
and every corresponding eigenfunction at depth n vanishes on V n -\ . 

The proof of this proposition can be found in pQ and in a large part depends 
on the Schur Complement formula and the eigenvector equations from previ- 
ously in the section. In this source it is shown how the eigenspace projection 
matrices are constructed, these matrices are related to the map — (D— z)~ 1 C and 
project eigenvectors on V n onto V n —\. This proposition is used constantly as we 
recursively extend eigenvalues from the V\ approximation, which are calculated 
explicitly, to the full fractal F m . A necessary condition for cr(M m „) capturing 
all of a (A m ) is that the spectral dimension is less than two [14] however it is 
known that for p.c.f. fractals this is always the case. 
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4 The Three-Branch Tree 



We analyze the spectrum of the Laplacian, A3, on the three-branch tree as a 
special case to demonstrate the calculations involved with spectral decimation 
as found in [TJ[2] while the approximating graphs and matrices are still relatively 
small. For this section assume that m = 3 throughout with will drop the m 
from the notation. Spectral decimation uses the spectra of Mo, Mi, D, and the 
Schur Complement of Mi — Iz yielding the spectrum of M„ as the output. We 
will then take the limit as n grows to infinity to find a (A3). The process will 
also be used in the next section as we consider the general case of the m-branch 
tree. 

First we consider Vq, which is simply a triangle (Figure [T]). The Laplacian 
matrix for Vq is 



M 



Since Vq is a complete graph there are no zeros in this matrix. 
The eigenvalues of Mq written with multiplicities are: 



<r(M ) 



3 3 
2' 2 



Next we consider the Laplacian on V± , which is a set of three triangles with 
a common vertex (Figure [T]) . There are seven vertices in V\ , and the depth-one 
Laplacian matrix is: 



Mi = 



V 





1 




1 

2 



1 

6 



\ 



1 / 



Recall from (|3.1j) that A will be the identity matrix I3 . 

The Schur Complement is S(z) — (A — z) — B(D — z)^ 1 C, as stated earlier. 
With these small matrices S(z) can easily be computed: 



S(z) 



( 

V 



(2.z-3)(6z -6z+l) 
6(z-l)(2z-l) 
2z-3 
12(z-l)(2z-l) 
2z-3 



12(z-l)(2z-l) 

From this, we have 



2z-3 
12(z-l)(2z-l) 
(2z-3)(6z 2 -6z+l) 
6(z-l)(2z-l) 
2z-3 
12(z-l)(2z-l) 



2z-3 
12(z-l)(2z-l) 

2z-3 
12(z-l)(2z-l) 
(2z-3)(6z 2 -6z + l) 
6(z-l)(2z-l) 



(z) = -(m - l)5i,2 



3 - 2z 



6(-l + z)(-l + 2z) 
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and 



R(z) = 1 



Si,! 

cj>{z) 



6z — 6z 



as defined in [U [2] . These two functions allow us to use Proposition 13.21 to 
determine the multiplicities of eigenvalues for M n , without writing the matrix 
M n down. These recursive formulae utilize the spectrum of M\ to recursively 
list all eigenvalues of M n . 

The eigenfunction extension map that fills in the interior values when eigen- 
vectors are extended is: 



/ 



-2+3z 



{D-z)- 1 C = 



3( 




-1+22) 


6( 


-l + 2)(-l + 22) 

-2+32 


6( 


-1+2) (-I+22) 


6( 


-1 + 2)^- 


-1+22) 


3( 


-l + 2)(-l + 22) 


6( 


-l + 2)(-l + 22) 

-2+32 


6( 




-1+22) 


6( 


-1 + 2)^-1 + 22) 


3( 


-1+2)^-1+22) 




3(-l+ 


22) 




3(-l+2z) 




3(-l+2 Z ) 




Figure 3: The graph of R(z) for the three-branch tree 

Now we list the eigenvalues of Mi, written with multiplicities, to provide a 
check for the calculations. 

, [3 3 3 3 1 1 1 

With corresponding eigenvectors of {-1, -1, -1, 0, 0, 0, 1}, {0, 0, -1, 0, 0, 1, 0}, 
{0, -1, 0, 0, 1, 0, 0}, {-1, 0, 0, 1, 0, 0, 0}, {-1, 0, 1, -1, 0, 1, 0}, {-1, 1, 0, -1, 1, 
0, 0}, and {1, 1, 1, 1, 1, 1, 1} respectively. 

To find the exceptional set we need the eigenvalues of D, which are written 
with multiplicities: 

a(D) = [1,1,1,1] , 

with corresponding eigenvectors of {-1, 1, -1, 1}, {-1, 0, 1, 0}, {-1, 1, 0, 0}, and 
{1, 1, 1, 1} respectively. 

The equation <p(z) — has only one solution, {§}, so the exceptional set is: 

E(M , Mr) = a(D) U {z : <j>{z) = 0} = \\, 1, \\ . 
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We can find the multiplicities of these exceptional values by using Proposi- 
tion [372j For the value of z = §, which is in cr(D), is not a pole of (j>(z), and 
<p(z) — 0, we use Proposition 13 . 21l5| ) to find the multiplicities: 

multo(j) = 2, 
multi(|) = 4, 
mult 2 (§) = 10. 

For the values of z = 4 and z = 1, which are poles of ^(2) and in <r(D), we 
use Proposition 13 . 2l| 3"j) : 

mult (|) = 0, multo(l) = 0, 
multi(|) = 2, multi(l) = 0, 
mult 2 (|)=4, mult 2 (l) = 0. 

For the value of z = 0, since ^ E(Mo, M), we use Proposition l3.21 fTj) to 
determine that mult„(0) = 1. 

Using the inverse function we can calculate the ancestors of z = 4 

to be i?-^ 1 / 2 ) = ^fi- From Proposition OP) , since 2^ £ E(Mo,M 1 ), 
their multiplicity at depth n will be the multiplicity of fl( 3± 6 ) = | at depth 
n-1. 

Table[T]shows the ancestor-offspring structure of the eigenvalues of the three- 
branch tree. The symbol * indicates an ancestor of 3± 6 , calculated by the 
inverse function R~ 1 (z) computed at the ancestor value z. By Proposition ^ . 2lf Tj) 
the ancestor and the offspring have the same multiplicity. The empty columns 
represent exceptional values. If they are eigenvalues of the appropriate M n , 
then their multiplicity is shown in the right hand part of the same row. 



z£a(M ) 
multo (z) 



1 


3 
2 

2 






2Gcr(Mi) 
multi (z) 



1 


1 


1 
2 

2 


3 
2 






zea{M 2 ) 
mult2 (2) 




1 


1 


3± 
( 

2 


2 


1 

2 


3 
2 




zea{M 3 ) 
mult3 (2) 




1 


1 


* 
2 


* 
2 


* 
2 


* 
2 


3±^6 
6 

4 4 


1 
2 


3 
2 

2S 



Table 1: Ancestor-offspring structure of the eigenvalues on the three-branch tree 
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Proposition 4.1. The dimension of the function space on V n , dim„, as was 
commented in Remark \3.1l is the same as the number of points in V n . For the 
three-branch tree, this is: 

dim n = 1 + 2-3". 
Likewise, mult„(z) is given as follows for n ^ and 1 ^ k ^ n: 

mult„(0) = 1, 
mult„(l) = 0, 
mult„(|) = 1 + 3", 
mult„(f) = l + 3 n " 1 , 
mult n (^5) = 1 + 3"~ 2 , 

mult n (i?- fc (|)) = 1 + 3"~ fc . 

Proof. For F$, the number of vertices begins with 3 in Vq then increases by 
3«-i _|_ 3" f or eacn subsequent level. Therefore, the partial sum of this sequence 
yields dim„ = 1 + 2-3". In Proposition ETJ we will prove dim„ for the general 
case. 

By applying Proposition 13 . 21 to the eigenvalues of Mi, and by using dim„ = 
1 + 2 - 3", 4>{z), and R{z), the multiplicities of a(M n ) are calculated inductively. 

□ 

Lemma 4.2. The spectrum of A3 is [J i>0 \R~ l (§)} U {0}. This spectrum is 
bounded, c(A3) C [0, 1) U {§}, and accumulates to the Julia set of R(z). 

Proof. From Figure [3] we can see that R(z) is a parabola with zeros {0, 1} and 
vertex (h, |). By applying inverses of R(z), we have R~ 1 (z) e [0, 1] so long as 
z e [0, §]. Since _R _1 (0) = {0, 1} and mult„(l) = 0, the only ancestor of will 
be itself. Since Ui> { R ~ l (I)} - (0> 1) u {§}- we have 

a(A 3 ) = U{i?-'d)}u{0}C[0,l)U{f}. 

The statement about the Julia set is a special case of the discussion in |5j 
Chapter 14. □ 

5 The m-Branch Tree 

This section focuses on the use of spectral decimation on the m-branch tree 
and its Laplacian, A m . We follow the same process as in Section [U however 
the inverse of (D — z) is more difficult to compute since it must be done for 
arbitrary m so we must use a different method for this computation. 

Proposition 5.1. For F m , the number of vertices in V m ^ n , which is equal to 
the dimension of the function space on V m ^ n (Remark \3.1)) . is 

dim OT ,„ = 1 + (m - l)m", 
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where m and n are integers with to ^ 3 and n ^ 0. 

Proof. Beginning with V m fi in F m , dim m _o = to. At each subsequent level, 
the number of vertices introduced by the defining iterated function system is 
(to- 2)m" +m n_1 . Thus, 



dim r , 



= m + ^](m-2)m , + m t 



i=l 



For any level n ^ 0, the partial sum of the sequence is dim„ = l + (m— 1)to™. □ 

Proposition 5.2. For £/ie depth-zero Laplacian matrix, mult (0) = 1 and 
mult (^ T )=m-l. 

Proof. For F m , we will have an to x to depth-zero Laplacian matrix of the form 



/ 1 — — 

m — 1 m — 1 



M 



m,0 



m— 1 

-1 -1 
m— 1 m— 1 



-1 
m— 1 
1 



^ m— 1 m— 1 m— 1 



m — 1 ' 

-1 
m — 1 

-1 
m— 1 



/ 



We are looking for the values of z and vectors, v, such that M mi0 w = zv. 
Since M m .o has a row-sum of zero v = {1, . . . , 1} G M m will be an eigenvector 
of M TO:0 with corresponding eigenvalue z = 0. So this is our first pair. 

Consider u = {0, . . . , 0, 1, —1, 0, . . . , 0} G R m , where 1 and -1 occur consec- 
utively and all other elements of v are 0, we get a corresponding eigenvalue 
z = -^-r. There are to — 1 of these vectors, so multnf— = m — 1. 

These vectors are linearly independent by standard methods, so they form 
a full set of eigenvectors for M m , - □ 

For the m-branch tree, M m .i is a square matrix with to 2 — to + 1 rows and 
columns. Take A to be the to x to identity matrix, then £> is a square matrix 
with (to — l) 2 rows and columns. Likewise, B will have dimension m x (m — l) 2 
and C will have dimension (to — l) 2 x m. Let 



5 = 



( m-1 








\ 



V o 



-1 



TO 



We will then have: 



M, 



m , 1 



,4 


B ~ 


C 





J 


(5 


m — 1 




J 










S 


6 




-l 


-l 


m — 1 
-1 l 


m(m— 1) 


m(m — 1) 


m(m— 1) > 
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In the last row and column the entries are row and column vectors with the 
indicated values in every entry. The lower right hand entry is a single entry 
however. 

Proposition 5.3. For the m-branch tree, cr(D) = ^rx, } with mul- 

tiplicities mult oiz^r) = m 2 — 3m+l, muhW— 2 TT ) =m— 1, and mult^f— rr) = 
1. 

Proof. We can observe that D has a row-sum of ^Uj- . Therefore, 

«i = {l,...,l} 

will be an eigenvector of D with corresponding eigenvalue Ai = ^zr[- 
We state the remaining eigenvectors and their eigenvalues. If 

«2a = {-l,...,-l,0,...,0,l}, 

where -1 occurs in the first m entries of the vector, v 2a will be an eigenvector of 
D with eigenvalue A 2 = ^rj- 
Consider 

v 2b — { x l> x 2i ■ ■ ■ -i x mi x li ■ ■ ■ i x m >0}i 

where Y^=i^ x j — 0, n °t e that the sum 1S across the upper index. If wc choose 
x\ = — 1, for some i > 1, x\ = 1, and all other = 0, v 2b will be an eigenvector 
of D with eigenvalue A 2 = — ^ . Since the 1 can occur in one of m entries, and 

& z m— 1 1 

the corresponding -1 can occur in one of m — 3 entries, v 2 & occurs in m(m — 3) = 
m 2 — 3m linearly independent variations, and consequently A has a multiplicity 
of m 2 — 3m + 1 (including the vector v 2a ). 
Consider 

V3 = {ci, C 2 , . . . , Cm, Cl, . . . , Cm, . . . ,0}, 

where c « — 0- If we choose c\ — —1, some a = 1, and all other Cj = 0, V3 

will be an eigenvector of D with eigenvalue A3 = -^zi ■ Since there are m — 1 of 
these vectors, A has a multiplicity of m — 1. 

It can be checked that v\, v 2a , V2b, and w 3 are linearly independent, so these 
vectors form a spanning set of eigenvectors. □ 

Proposition 5.4. The two entries of the Schur complement used in calculating 
R(z) and 4>{z) are: 

&,! = (!-*)- ™--(--l) 2 



to(1-z(to-1))(2-z(to-1)) 
and 

m — z(m — 1) 

Ji,2 ~~ 



m(m - 1)(1 - z(m - 1))(2 - z(m - 1)) 
/or any m ^ 3. 
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Proof. To find the Schur Complement, one must invert the increasingly large 
matrix D — z, where D is an (to— l) 2 x (to— l) 2 matrix. This is most conveniently 
done by forming a spectral resolution of D. Since D has 3 distinct eigenvalues, 
we can express D as: 

D = Ai • Pi + A 2 ■ P 2 + A 3 • P 3 , 

where Ai is an eigenvalue of I? and Pi is the projection matrix onto the eigenspace, 
Ei, corresponding to the eigenvector A^. This can be done since the eigenvectors 
are all linearly independent [13] . Once the projectors are written, then, 



(D-zy 



1 



Ai — z 



P 



1 



A, 



P2 + 



1 



A, 



P 3 



To find the projection matrices, we must find matrices that have the following 
properties: 

1. map each eigenspace to itself, Pi(Ei) = Ei; 

2. Pi = ft; 

3. Pi ■ v — when v £ © Range(Pj); and 

4. ft + P 2 + P 3 = I. 

First we consider Pj and P3, both of which are (to — l) 2 x (m — l) 2 matrices: 



1 



ft 



1 



to(to — 1) 



1 TO 



1 TO 



(m— l) 2 X (m — l) 2 



1 



to(to — 2) 



J • 

V • 



j \ 

j 
00/ 



(m-l) 2 x(m-l) 2 



where J is the mxm matrix 



/ to - 1 -1 
-1 to-1 



J = 



V -1 



-1 



-1 \ 

-1 



1 / 

' mxr 



Then we set P2 = I — (Pi + P3) to satisfy the fourth requirement for the 
projection matrices listed above. While adding Pi and P3, we obtain a new 
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to x to matrix, K. Subtracting Pi + P 3 from I, we have: 



(I-K -K 

-K I-K 



P2 



K 

-l 



-K 



-K 
-K 



m— 1 



where 



X = 



-1) m(m — 1) 



/ m 2 — m — 1 
-1 



-l 

m— 1 
-1 m-2 



7-X 

ro-2 / 

*("»-!) m-1 ' (m-l) 2 x(m-l) 2 



m(m — l)(m — 2) 



V 



-1 

to 2 — m — 1 



-1 



-1 \ 

-1 

m 2 - m - 1 J 



It is not difficult to show that these are projecting onto the correct subspaces 
by checking their action on the eigenvectors. 

Now that we have our projection matrices and their corresponding eigen- 
values, and since matrix multiplication is distributive, we can write the Schur 

3 

Complement as: S(z) = (A — z) — £ -y^—^BPiC. The 1, 1 and 1, 2 entries are 

i=i 1 2 

the only ones we need to calculate <j>(z) and R(z). The relevant entries from the 
summands are: 



T^-BP lCl ., 

J^BPid.2 



Pi 



Pi 



Pi 



o 



(l-z(m-l)) 

1 



i(l-z(m-l)) 
3 

I 

i=l 



m(2-z(m-l)) 
-(m-2) 



m(m-l)(2-z(m-l)) 

Since (B(£> - z)" 1 ^!,! = £ ^(BPjCOi.i and (Bp - ^)- 1 C) 1 , 2 



£ T J_(BP j C , )i, 2) we have 
i=i 



and 



S M = (l-z)- 



m — z(m — l) 2 



to(1-z(to-1))(2-z(to-1)) 
to — z(m — 1) 



m(m - 1)(1 - z(m - 1))(2 - z(m - 1)) ' 



□ 



Corollary 5.5. For any m ^ 3, ifte function 4>{z) = — (m — l)Si,2 

m — (m — l)z 



m(2 — (to — l)z)(l — (to — l)z) 
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Figure 4: The graph of R{z) for the m-branch tree 

Corollary 5.6. For any m 3, the function R(z) = 1 — is: 

R(z) — 2mz — m(m — l)z 2 . 

Theorem 5.1. for any m ^ 3, n ^ 0, and 1 ^ fc ^ 77, i/ie depth-n Laplacian 
matrix will have eigenvalues such that: 

mult n (0) = 1, 

mult n (^_) = 1 + (to - 2) • to™, 
mult„(fl- fc (^ I )) = 1 + (m - 2) ■ m"- fe . 

In fte {imtt <r(A m ) = U t > {^(t^t)} U {0} C [0, ^) U {^}, 

where the inverse images of n "^ 1 accumulate to the Julia set of R(z). 

Proof. We first define the exceptional set, E(M m ^, M m ) = a{D) U {z : <j>{z) = 
0}, as in Section [3] Since <fi(z) = only when z = ^zj, we have 

f m 2 1 1 
E(M mfi , M m ) = -, -, \ . 

( 777 — 1 777 — 1 777 — 1 J 

We can find the multiplicities of these exceptional values using Proposition [3T2J 
For z — zprzr ; which is in a(D), is not a pole of 4>(z), and 0(z) = 0, we use 
Proposition 13. 2t [5|) to find its multiplicities: 

multnC^j) = 1 + (777 - 2) • 777". 

For z = rr^r, which is in cr(D), is a pole of an d R'( z ) 0; we use 

Proposition 13. 2t f3"|) to find its multiplicities: 

mult n (^) = 0. 

For z = which is in a(D), is a pole of <fr(z), and = 0, we use 

Proposition 13. 2l[ 6|) to find its multiplicities. Note that multo( ) = so: 

mult„(^_) = multn-i^) = 1 + (777 - 2) • m"" 1 . 
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When z = 0, since z £ E(M m ^ : M m ), we use Proposition 13. 2p | to find its 
multiplicities. Given multo(O) = 1 and R(0) = 0, 

mult„(0) = 1. 

Since i? _1 (0) = {0, m _i } and mult n ( _-, ) = 0, there are no additional 
ancestors from 0. However, as n increases, one must consider R~ k (^hj) for 
1 < k < n. Since R~ k (^j) £ E(M , M) when k > 1, we use Proposition l3T2"lTTj) 
to determine that 

mult„ (R- k (yS=x)) = mult n _ fc (^) = 1 + (m - 2) • m"" fc . 

As n — » oo, the only ancestor from z = is itself. However, since R{z) has 
zeros at z = {0,^-} and maximum at R(^j) = R^ 1 (z) e [0, 7 ^ rT ] 

so long as z e [0, ^j], as seen in Figured] Since [j i>0 (^t) } - 

(0,^ T )U{^ I },wehave 

a(A m ) = |J {i?-(^r)} U {0} C [0, ^) U {^}. 

As before, the inverse images accumulate to the Julia set of R(z) as discussed 
in [5]. □ 

Table [2] shows the ancestor-offspring structure of the eigenvalues of the 
m-branch tree for the first few Laplacian matrices. In the table, we denote 
R~ 1 (-^i) = ip and R^ 1 (ip) = *. Blank entries denote exceptional values and 
if they have non-zero multiplicity they are added at the right hand end of the 
table. 

Proposition 5.7. The eigenvalues and multiplicities established in Theorem 
\5.1\ revresent all the eigenvalues of M m n . That is, at any level n, the number 
of eigenvalues is equal to dim n . 

Proof. We can conclude from Theorem 15.11 and Table [2] that all eigenvalues for 
any level n in F m lie in either [0, ^zj) or {— The number of eigenvalues 
in the depth-n Laplacian of the m-branch tree will be 

n 

mult„(0) + mult„(^ T ) + mult* (V 1 ^)) 

n 

= 1 + (1 + (m - 2)m n ) + j^(2 <_1 (l + (m - 2)™""*)) 

i=l 

= 1 + (1 + (m - 2)m n ) + (2™ - 1) + (m n - 2") 
= 1 + (to — 1)to™ = dim„. 



The equation holds because R~ 1 (z) will yield exactly two values for all z € 
(0, rr^r)- Since dim n = 1 + (to — 1)to™, we have found every eigenvalue for the 
depth-n Laplacian matrix. □ 
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Z£i(T(JV1q ) 

multo(^) 



1 


m 
m — 1 

m — 1 










multi (z) 



1 


2 

m-1 




m 

m 


L 

-1 
- 1 




m 
m— 1 

(m-1) 2 






zea(M 2 ) 
mult2 (z) 




1 


2 

m — 1 




m — 1 


m 


2 

- 1 




m 

(m - 


-1 

"If 


in 
m— 1 

tti 3 - 2m 2 + 1 




zEa(M 3 ) 

mult3(2) 




1 


2 

,71-1 






* 

m — 1 


* 

m - 1 


* 

m - 1 


* 

m-1 


in 

(m - l) 2 


^2 

(m - l) 2 


1 

m— 1 
m 3 - 2m 2 + 1 


m 
m — 1 

™ 4 - 2m 3 + 1 



Table 2: Ancestor-offspring structure of the eigenvalues on the m-branch tree 



It is of some interest as well to have some idea of the geometric properties of 
the fractal as well as its vibrational behavior. To this end we need more detailed 
information about the generating iterated function system than the qualitative 
description given in the original description of the fractals. 

Proposition 5.8. For any n ^ and m ^ 3, the common contraction factor 
c m for all functions in the IFS as measured in the effective resistance metric is: 

_ 1 

Cm - -. 

Proof. The idea behind this is that if two elements of Vo are separated by 
resistance one, then those same points viewed in V\ still should have resistance 
one between them. So you pare off the uninvolved branches and get that they 
are now at opposite ends of a chain of two copies of Vq each of which would 
have equal resistance by symmetry so they must have resistance one half. So 
the distance between neighboring points goes down by a factor of one half each 
time the level of approximation is increased. For a more formal treatment of 
these types of calculations see [SJ [TJ] . □ 

Once the contraction factor is known it is a direct corollary to calculate 
the dimension of a given to— branch tree. These fractals satisfy the open set 
condition which is a prerequisite of the formula in this proposition. 

Proposition 5.9. Since the m-Branch Trees satisfy the open set condition and 
are generated by to contraction mappings the Housdorff dimension s m is given 
by: 



—2 log(TO) 
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Proof. A formula for the Hausdorff dimension of a fractal generated by an IFS 
with a finite number of functions is given in [5]. The formula is: 

m 
i=l 

In our case Cj = h so the calculations are simple and yield 

—2 log(m) 
Sm = i 7T~\ • 

□ 

This confirms the reasonable intuition that as the number of branches grows 
so does the dimension of the space. 



6 Towards the Infinite-Branch Tree 



With the full spectrum of the m-branch tree established in Theorem 1 5. 11 we can 
observe what happens as the number of branches goes towards infinity. From 
Propositions 13.21 and 15.11 and Theorem 15. 1[ it is clear that the multiplicities 
of eigenvalues and dim mj „ both go to infinity. However, we can observe what 
occurs to the density of the spectra as m — > oo. 

As a tool to describe what happens to the spectra of the A m as to grows 
we introduce a measure on the complex plane supported on er(A TO ). For a 

z € Rr k (^^zjj then for some n large enough mult„(z) is positive. Also is an 

eigenvalue with mult„(0) = 1 for all n. Points, z, of these two forms are dense 
in a(A m ) so we define the measure as follows 



hm — 

n— >oo dirrin * — ' 

n z£a(A m ) 



mult n (z)<5 z 



Now we consider the weak limit of these measures and the support of the limit 
measure. 

Theorem 6.1. As limm^oo <r(A m ) = {0,1} in the Hausdorff metric. Also 
Iim m _ >oo K m — where the limit is in the weak topology. 



Proof. From Theorem 15.11 we have that 



u * 



i>0 



m — 1 



c 



0, 



m — 1 



U 



TO — 1 



and 



a(A m )=\j{R-\ 1 ^ I )}u{0}C = 

i>0 



{nl-l} 
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It is easy to see that the only points in the complex plane about any open ball 
will intersect infinitely many of the cr(A m ) are the points {0} and {1}. 

Since mult„(— ^) = 1 + (m — 2)m" and dim„ = 1 + (m — l)m n , the density 
5F of in the spectrum of M n as m — > oo is 

(to \ 1 + (m — 2)m™ to™ +1 — 2to" + 1 ^ 
to- 1 J ~~ 1 + (to - 1)to" ~~ m «+i - m « + l ^ m ^°° 

And similarly the density of all other eigenvalues which are the ones tending 
towards zero is O(^) as to — > oo. 

□ 
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